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Abstract 

We study the low energy quantum spectra of two-dimensional rectangular bil- 
liards with a small but finite-size scatterer inside. We start by examining the 
spectral properties of billiards with a single pointlike scatterer. The problem 
is formulated in terms of self-adjoint extension theory of functional analysis. 
The condition for the appearance of so-called wave chaos is clarified. We then 
relate the pointlike scatterer to a finite-size scatterer through the appropriate 
truncation of basis. We show that the signature of wave chaos in low energy 
states is most prominent when the scatterer is weakly attractive. As an illus- 
tration, numerical results of a rectangular billiard with a small rectangular 

scatterer inside are exhibited. 
5.45.+b, 3.65,-w 
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I. INTRODUCTION 



The two-dimensional billiard is an appropriate tool for examining the generic features 
of dynamical systems because of the wide range of dynamical behaviors going from the 
most regular (integrable) to the most irregular (chaotic) depending on the geometry of its 
boundary Although no mathematical proof exists, it is widely believed that fingerprints of 
the regular or irregular nature of the classical motion can be found in statistical properties 
of the quantum spectrum both in energy levels and wavefunctions. The integrable systems 
such as the circular, elliptic and rectangular billiards obey the Poisson statistics [T]|J, while 
the predictions of the Gaussian orthogonal ensembles describe the chaotic systems such as 
Sinai's billiard and Bunimovich's stadium [[HQ. 

Besides these extreme classes, there is an intermediate category called pseudointegrable 
or almost integrable system. The typical examples are staircase billiards || , rational 
polygons [|^,|9|-|TT|j and the billiards with pointlike scatterers inside an integrable one fi"2|-|T7 . 



The nature of classical motion in pseudointegrable systems can be considered as integrable 
in the sense that unstable trajectories are of measure zero. However, several numerical 
studies have revealed that under a certain condition, the quantization induces the chaotic 
energy spectra, which can be regarded as a counterexample for the correspondence between 
the energy spectra and the underlying classical motion. This phenomenon is named wave 
chaos because its origin is the wave-like nature of quantum motion. For the billiards with 
pointlike scatterers, this phenomenon has been understood in terms of the quantum breaking 



of the classical scale invariance |17| . However, the notion of the pointlike scatterer is a 
mathematical abstraction whose relevance to the physical system is far from straightforward. 
It is therefore highly desirable to show that the nature of wave chaos remains intact for the 
case of finite-size scatterer. This is the prime motivation of the work we report in this paper. 
We restrict ourselves to considering the billiards with a small but finite-size scatterer inside 
an integrable one. It is shown that, even if the scatterer has a substantial size, the signature 
of wave chaos are observable under a certain condition. The condition is consistent with 
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a previous numerical observation in Sinai's billiard, which indicates that there exists no 
signature of wave chaos when the scatterer is repulsive. 

The paper is organized as follows. We start by examining the spectral properties of 
pseudointegrable billiards with a pointlike scatterer in Sees. II and III. In the former section, 
we give an accurate description of a pointlike scatterer, based on the self-adjoint extension 
theory in functional analysis. A special emphasis is laid on the distinction between symmetry 
and self-adjointness. We give an intuitive explanation of the necessity of the self-adjointness 
for quantum mechanical Hamiltonians. Applying a general prescription for extending a 
symmetric operator to a self-adjoint one, we derive the Green's function for the billiards with 
a single pointlike scatterer. A renormalized coupling constant of the scatterer is defined in a 
natural manner, although its physical meaning is somewhat unclear at this stage. Based on 
the formulation in Sec. II, we discuss the quantum spectra of the billiards with a pointlike 
scatterer in Sec. III. The general condition for the appearance of wave chaos is clarified. In 
Sec. IV, we consider a Dirac's delta function potential with a truncated basis. Although it is 
not well-defined with a full basis in two-dimensional billiards, we can deduce certain physical 
contents from the truncated system, which serve as the basis for the ensuing discussions. 
The findings in Sees. Ill and IV are applied to the cases of a small but finite-size scatterer in 
Sec.V. It is shown that the wave chaos is manifest at low energy with attractive potentials. 
As an illustration, we give a numerical example for a rectangular billiard with a small 
rectangular scatterer inside. We also show a proper procedure for the zero-size limit of the 
finite-size scatterer, which is consistent with the self-adjointness of the Hamiltonian. This 
clarifies the physical meaning of a renormalized coupling constant defined in Sec. II. We give 
the conclusions in Sec. VI. 
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II. QUANTUM MECHANICAL FORMULATION OF PSEUDOINTEGRABLE 
BILLIARDS WITH A POINTLIKE SCATTERER 



In spite of the apparent simplicity, careful treatments are required for the definition of 
quantum billiards with a pointlike scatterer inside. Several methods are known for this 



purpose [18]. Here, we adopt the one based on the self-adjoint extension theory in func- 



tional analysis |B5|. We briefly summarize the formalism stressing the necessary points for 
the discussions in following sections. We also mention the physical reasons why quantum 
mechanics requires the self-adjointness for Hamiltonians. 

Let us consider an integrable billiard of area S with the Dirichlet boundary condition 
such that wavefunctions vanish on the boundary. The eigenvalues and corresponding eigen- 
functions are determined by the stationary Schrodinger equation; 

H (f n (x) = -^A(p n (x) = E n if n (x){n = 1, 2, 3, • • •), (1) 

where x is the coordinate vector in two-dimensional space, and M is the mass of a particle 
moving in the billiard. We adopt natural units throughout this paper; h = l,c = 1. This 
leaves a single independent unit among mass, energy and length. The domain of Hq is 
D(Hq) = H 2 (S) fl Hq(S) in terms of the Sobolev spaces. (We denote the domain of the 
billiard by the same symbol as the area unless there is danger of confusion.) In the following, 
we assume that the unperturbed billiard has no degeneracy. The Green's function of this 
system is given by 
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C<°>R^) = £^%^ (2) 

71=1 Z n 

where z is an energy variable. Suppose that a pointlike scatterer is placed at x = x$ inside the 
billiard. Following the self-adjoint extension theory in functional analysis, the eigenvalues 
and the corresponding eigenf unctions of the perturbed system are calculated as follows. 

The first step for this purpose is to remove the relevant scattering point Xq by restricting 
Hq to T = Hq\D(T) with the domain 
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D(T) = Mx) E D(H ) | <p(x ) = 0}. (3) 

By using integration by parts, it is easy to see that T is symmetric (Hermitian). However, 
the domain of T*, the adjoint operator of T, is not identical to D(T) and indeed larger than 
D(T); 

D(T*) = D(T) © Ran{T - X) L © Ran(T - A) x , (4) 

where the orthogonal complement to the range of T — A corresponds to the deficiency sub- 
space of T given by fl20j 



Ran{T- X) L = {<p(x) \ tp{x) = cG {0) {x, x , A), c G C} (5) 

with an arbitrary complex number A (ImX ^ 0). It follows from D(T*) ^ D{T) that T is 
not self-adjoint. A symmetric operator is self-adjoint if and only if D(T*) = D(T). 

The two following facts serve to understand the reason for imposing self-adjointness on 



quantum mechanical Hamiltonians [21,22|. One of them is based on the assertion 



Ker{T* - A) = Ran{T - A) . (6) 



This means that T* has a complex eigenvalue A with eigenspace Ran(T — A) , namely, the 
closure of Ran(T — A) -1 . If and only if T is self-adjoint, Ran(T — A) -1 = 0, i.e., Ran(T — A) 
is the entire Hilbert space X = L 2 (S), the set of all square- integrable functions over S. This 
is one of the reasons why quantum mechanics, which does not allow complex eigenvalues 
for observables, requires self-adjoint Hamiltonians. The second fact is closely related to the 
unitarity of time-evolution operators. The time-evolution operator U (t) is constructed from 
T as follows; 

U(t) = e" m = lim (1 + —)~ n = lim {(-)-"(T - -)""}. (7) 

It is easy to see that, for any symmetric T, an operator (T — A) -1 with pure imaginary A ^ 
is defined from Ran(T — A) to D(T) with || (T — A)~V|| < ||<^||/|/mA| for any (p e Ran(T — A) . 
If Ran(T — A) is not all of X, then in general, (T — A) _ra can be defined on smaller and 



smaller spaces as n increases. This leads us to a natural condition Ran(T — A) = X, 
implying that time-evolution operators which preserve probability can be constructed only 
from self-adjoint operators. 

The above considerations indicate that we have to extend D(T), or equivalently, restrict 
D(T*), in order to construct a self-adjoint operator from T = H \D(T). This is indeed 
possible since dimRan(T — A) -1 = dimRan(T — A) -1 . Because both the deficiency subspaces 
are one-dimensional, a one-parameter family of self-adjoint extensions of T exist. According 
to a general prescription for extending a symmetric operator to self-adjoint ones [23|, All 



the self-adjoint extensions of T are given by Hg = Ho (0 < 9 < 2ir) with the domain 

D(H g ) = {ip(x) | ij){x) = (p{x) + cG {0) (x,x ; A) - ce ie G^(x, x ; A); 

<p(x) G D(T), G (0) (x, x ; A) e Ran(T - A) x , G (0) (x, x ; A) G Ran(T - A) x , 
cgC}. (8) 

Eq.(pp with Eq.(^|) indicates that the operation of Hg on ip{x) is given by 

(H e i))(x) = (H <f)(x) + XcG {0 \x, x Q ; A) - Xce id G {0 \x, x ; A). (9) 

The operator Hg is regarded as the Hamiltonian of the perturbed system with a pointlike 
scatterer. Notice that, although we have Hg = H on D(T), D(Hg) is substantially larger 
than D(T); D(T) C D(Hg) C D(T*). We recognize from the second and third terms of 
ip{x) in Eq.(P) that the appropriate boundary condition around the scatterer is specified 
after the extension. It will be shown later that the value of 9 is related to the strength of 
the pointlike scatterer. 

With the aid of the resolvent equation, we can write down an equation for the Green's 
function of the system with the Hamiltonian Hg as 

G e (x, y; z) = G (0) (x, y; z) + G {0 \x, x ; z)T e (z)G^ (f , V\ z), (10) 

where the transition matrix (T-matrix) Tg(z) is calculated through 
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1 ~ eie 

9 ^ ~ (z-X)J (x, x o; (£, f o; A)dx - e* e (z - A) / G(°) (f, x Q ; z)G® (x, x ; \)dx 

(11) 

After substituting Eq.@ into Eq. flTT|) , and with a somewhat lengthy but straightforward 
calculation, we have 

T e (z) = (v e 1 -G(z))-\ (12) 



where 



-1 _ IM Sin (l + ar ^ A ) ^ ^n(^o) 2 , lq x 
sm 2 n=ll- C 'n -A l 



oo IF 

G(s) = £ ^ B (* )*(__ + ttt-V). (14) 

n= l 2 — £/ n |£/ n — A | 

It follows from Eq.(|12|) that the perturbed eigenvalues are determined by 

G{z) = vj\ (15) 

The solutions of this equation, z n , (n = 1, 2, 3, • • •), correspond to the poles of the T-matrix. 
The corresponding eigenfunction is the residue of the Green's function of Eq. fllOD at the 
pole; 

ij) n {x) =N n G^(x,x ;z n ) (16) 
where the normalization factor is determined by 

Although the T-matrix seemingly has two independent parameters, A and 9, this is not 
the case. This follows from the fact that |A| is regarded as a scale of mass in case of A being 
pure imaginary. Indeed, taking |A| as a unit of energy, we can fix A = +i, which makes all the 
relevant physical quantities to be dimensionless. Such treatment is justified by the assertion 
that a symmetric operator T is self-adjoint if the condition Ran(T — A) = Ran(T — A) = X 



is valid for some fixed A (JmA 7^ 0). This means that if and only if Ran(T ± i) = X, a 
symmetric operator T is self-adjoint. Thus, without any loss of generality, we can make the 
replacement of |A| — > 1, z — > \X\z, E n — > \X\E n and <p n — > \\\ip n in Eqs.([L3]), ( pT4| ) and (|15|). 



We then have 



= v \ (18) 



where 



_ sin9 ^ y?„(£ ) 2 



E^fV, (19) 



1 -cos# n=1 ~„ 

co IF 

G(s) = £ ^(f ) 2 (--r + -^t). (20) 

n=l 2 _ & n + J- 

We can regard as a coupling constant of a pointlike scatterer. It ranges over all real 
numbers as < 9 < 2n. 

In a previous publication ||16|| , we have referred to the coupling constant Vg as the bare 
coupling constant, denoted by vb- However, this is somewhat confusing, because vq should 
be considered to arise in a renormalization process for treating short-range singularities in 
a proper manner, which we often encounter in field theory. It might be more appropriate to 
call vg the renormalized coupling constant. 

Because the average level density of two-dimensional billiards is independent of energy z, 
each term in the parenthesis of Eq.(^) diverges when summed separately. The divergence 
disappears when summed together. This means that the second term of Eq.(|20|) plays an 
essential role to make the T-matrix well-defined. It is also noteworthy that the energy 
dependence appears only in the first term of Eq. fl2"0|) . This ensures the orthogonality of the 
perturbed eigenfunctions {ifj n (x)}, (n = 1, 2, 3, ■ ■ •). Namely, for m 7^ n, we have 



^ m (x)Mx)dx = N m N n G{Zn) _ G{Zm) = 0. (21) 



The orthonormal relations of the unperturbed eigenfunctions {(p n (x)}, {n = 1, 2, 3, • • •), are 
used in the first equality, and the second equality results from of Eq . (|18D . The authors of 
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Ref. [24 have discussed a rectangular billiard with a single pointlike scatterer and deduced 



an eigenvalue equation, Eq.(23) in Ref. |24|, which is seemingly analogous to Eq.(|18D with 
Eqs.([l9|) and (^) in this paper. However, their approach allows the coupling strength vg to 
be varied depending on the energy, and violates the orthogonality, Eq. (pip . 

III. QUANTUM SPECTRUM OF PSEUDOINTEGRABLE BILLIARDS WITH A 

POINTLIKE SCATTERER 

Equipped with the formulation in the previous section, we examine the spectral properties 
of pseudointegrable billiards with a pointlike scatterer. A special emphasis is placed on the 
condition for the appearance of wave chaos. 

An important fact is that each perturbed eigenvalue is isolated between two unperturbed 
ones. This follows from the fact that G(z) of Eq. fl20|) is a monotonously decreasing function 
of z on the interval between any two successive unperturbed eigenvalues, (E n , E n+1 ), covering 
the entire value between (— oo, oo). It is also clear that G(z) has a single inflection point 



on (E n , E n+ \). We have shown in Ref. [IS] that the disturbance by the pointlike scatterer is 
restricted to the eigenstates with an eigenvalue around which G(z) has an inflection point 
(See Figs. 2 and 3 in Ref. |16||). This is understood also from the eigenfunction, Eq.(|l6|) 
with Eq.(0); If a perturbed eigenvalue is close to an unperturbed one, the corresponding 
perturbed eigenfunction is not substantially different from the corresponding unperturbed 
one. 



In Ref. |16j , the condition for the inflection points has been deduced in a somewhat com- 
plicated manner, by making a truncation of the unperturbed basis according to the energy 
under consideration. Here, we refine the argument without introducing any truncation of 
the basis. 

Each inflection point of G(z), say z m , (m = 1, 2, 3, • ■ ■), is expected to appear, in average, 
around the midpoint on (E m , E m+ i); z m ~ (E m + E m+ i)/2. The contributions from < 
E n < z m and z m < E n < 2z m ~ E 2m on the summation of the first term in Eq.(]2Tj) are 
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canceled with high degree of accuracy. This allows us to estimate G(z) at the inflection 
point z m as follows; 



2 in 
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n=l "^n + 1 n=2m+l 



/•B2m E r°° 1 E 

Vo £ 2 + l JE 2m y ~z m -E E 2 + V 1 
~alnz m , (22) 

where we define a = p a v((Pn(xo) 2 ) in terms of the average level density of the unperturbed 
system, p av . The symbol (•) means an average among various n. Because p av = MS/(27r) 
according to the Weyl's theorem and ((p n (x ) 2 ) ~ 1/5* for a generic choice of x , we obtain 

« = Pav(Vn(x ) 2 ) ~ — . (23) 

It follows from Eqs.(|l8D and (p2|) that the effects of the pointlike scatterer on the quantum 
spectrum are observed in the eigenstates with an eigenvalue z such that 

alnz ~ Vg . (24) 

The condition of Eq.(p4[) corresponds to Eq.(60) in Ref. It is also noteworthy that a s- 
wave phase shift of the two-dimensional scattering problem with a single pointlike scatterer 
shows a similar logarithmic dependence on energy z Hl8| , ^5| . The origin of such energy 



dependence is closely related to the energy scale |A| introduced by the self-adjoint extension 
of a symmetric operator; Any function which depends on |A| has an energy dependence so as 
to guarantee the argument to be dimensionless. Although the nature of a particle motion in 
billiards is independent of the energy in classical physics, the singularity of the interaction 
induces an energy dependence of the observables after quantization. This can be considered 
as a typical example of quantum mechanical breaking of scale invariance, or scale anomaly 

Let us proceed to make an estimate of the width of a strip along the logarithmic curve 
of Eq.(^) on which the disturbance of the scatterer is observable. Using the derivative of 
G(z) at an inflection point z m , we can estimate the width, say A, as follows; 
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We have implicitly assumed in Eq.(^) that the unperturbed eigenvalues are distributed 
within a mean interval p~l in the whole energy region including negative energies. This 
assumption is quite satisfactory in this case, because the denominator of G'(z) is of the 
order of (z — E n ) 2 , indicating that the summation in Eq.(^) converges rapidly. With the 
aid of a relation 
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£ (2- - I) 2 " 8 ' (26) 



we arrive at 



G'(z m )\p-^7r 2 a, (27) 



which means 



A ~ a. (28) 

The estimate of Eq . (p8|) is justified from another perspective later in this paper. 

The above considerations lead us to the condition for the appearance of wave chaos. The 
effect of a pointlike scatterer with strength vg appears mainly in the eigenstates with an 
eigenvalue z satisfying 



i -— l 



cdnz|<A, (29) 



In other words, wave chaos occurs in the energy region which satisfies 

exp{— - 1}< z <exp{— + 1}. (30) 

OLVg ~ ~ OlVq 

The condition of Eq . (^9|) contains all the essential physics of wave chaos. From this, we can 
draw the following conclusions. 

(1) For any positive vg, wave chaos appears at the energy which satisfies Eq.(^), while it 
is hardly seen in other energies. If vg is larger than A , wave chaos appears only around the 
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ground state region. As the energy increases, it tends to disappear. On the other hand, if vg 
is substantially smaller than A -1 , wave chaos appears in the higher energy region specified 
by Eq.(|3(]). In the limit of vg — > +0, the system restores the integrability by pushing up the 
chaotic region to the infinite energy. (In this limit, the eigenvalue of the single eigenstate 
with an eigenvalue smaller than Ei diverges to — oo.) 

(2) Wave chaos does not appear at any energy if vg is negative. In particular, the system 
converges to the unperturbed one as vg — > — 0. 

(3) For any vg, the system behaves as integrable in the high energy limit, which agrees 
with our intuition that a pointlike scatterer has no effects on a particle motion in the classical 
limit. In other words, the quantum billiards with a pointlike scatterer have the nature of 
asymptotic freedom, which has been first discovered in the non-Abelian Gauge field theories. 

IV. FORMAL DESCRIPTION OF A POINTLIKE SCATTERER IN TERMS OF A 

DIRAC'S DELTA FUNCTION 

We have revealed various features of the two-dimensional billiards with a single pointlike 
scatterer from a general perspective in the previous section. ^From a practical point of view, 
however, we still need to identify the physical meaning of the coupling strength vg, which 
does not have a direct relation to physical observables. As a first step for this purpose, we 
show in this section that vg can be related to the physical strength of a pointlike scatterer 
through a truncation of the unperturbed basis. 

Let us consider a formal Hamiltonian with a Dirac's delta function potential 

H = H + v5(x-x ) (31) 

with a truncated unperturbed basis {{p n (x)}, (n = 1, 2, 3, • • • , N). We can regard v as 
the physical coupling constant of the pointlike scatterer. Notice that the Hamiltonian as 
Eq.fl3~T|) cannot be defined in the billiard problems with higher dimension than one because 
one cannot take the limit of N — > oo. For a moment, however, we will go further along with 
fixed N and return to this point later in this section. 
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The characteristic equation of H is given by 

det(zl - H - V) = 0, (32) 
where I is a iV-dimensional unit matrix. The matrix elements of Ho and V are given by 

(Ho)rrm = E n 6 mn (33) 

and 

(V) mn = v<p m (x )(p n (x ) (34) 

with m,n = 1, 2, 3, • • • , N, respectively. Using the linearity of a determinant and noticing 
that any two column vectors of V are linearly dependent, we obtain 

N N N 

det(zI-tt -V)=Y[(z-E n )-vY,Wl(xo) li(z-E m )}, (35) 

n=l n=l m=l 

where n' signifies that the term of m = n is removed from the product. By inserting 
Eq.fl3"5]) into Eq.(^), the characteristic equation can be rewritten as 

G{z) = v-\ (36) 

where 

n=l z ^ n 

For any v, Eq.(^) has non-degenerate solutions z n , (n = 1, 2, 3, • • • , N), giving a complete 
set of eigenvalues of the Hamiltonian Eq.(|3~l~D with the truncated basis. The corresponding 
eigenfunction is given by 

M*)=N„<t ^f°> (38) 

k=l ^ n fe 

with the normalization factor determined by 

k=1 [Zn - &k) 
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The eigenfunction, Eq.(^) has a form analogous to Eq. (|T6|) . For large N, the appropriate 
boundary condition around the scatterer is taken into account in a matrix diagonalization 
with a Dirac's delta function potential. 

In an analysis similar to the previous section, we recognize that the inflection points of 
G(z) are located around G(z) ~ at energy z ~ -Ejv/2, implying the appearance of wave 
chaos under the condition of v^ 1 ~ 0. We also expect that it disappears if |t> _1 |>A. In case 
of \v~ 1 \ ~ A, the average ratio between the diagonal matrix element of the delta function 
potential and the mean level spacing is estimated as 



<(V)nn> 



Pav 



| - 1. (40) 



This explains the reason why the logarithmic curve of Eq.(^) has a strip with width A; 
In case of |t> _1 |>A, the physical strength v is too weak to cause the mixing among a large 
number of unperturbed eigenstates. In this sense, we can regard If" 1 ! as a measure of 
distance to the wave-chaotic curve. We stress that | 1 1 ~ corresponds to the condition 
for the appearance of wave chaos in the energy region around z ~ -Eat/2. Both at the lower 
and higher energies, the value of G(z) at the inflection points is substantially different from 
zero, since the contributions on G(z) from the terms with E n < z and those with E n > z do 
not cancel. This fact gives rise to an interesting physical result in the following section. 

Eq.(p7|) allows us to relate the physical strength v to the formal strength vq introduced 
in Sec. II. In a similar manner as in Eq.(^), we get for z < En, 

N rp oo IF 

G(z)=G(z) + ^^(x o r-^- + E ^ f o) 2 (-^r + ^rr) 

n=l K + 1 n =N+l Z-E n E n + l 

rE N e r°° 1 E 

= G < 2 > + a U WTT dE + L ( 7^E + WTT )dE} 

~G(z) + aki(E N -z). (41) 
^From Eqs.([T8|), (|36|) and (|4lD , we obtain the relation between both the coupling constants; 



v. 



v' 1 + a\n(E N - z). (42) 



We note that the summation in Eq.(^) diverges as iV increases for any z ^ E n . This 
follows from the fact that a defined in Eq. (p3|) is independent of energy z. This underscores 
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the fact that the Dirac's delta function potential with a full unperturbed basis cannot be 
defined in two-dimensional billiard problems, contrary to in one dimension. In the latter 



to make the matrix diagonalization with strength v identical to the formulation with vq 
discussed in Sec. II. Thus it is possible to construct quantum mechanics of one-dimensional 
billiards with a pointlike scatterer directly in terms of a Dirac's delta function without 
relying on a somewhat complicated mathematical framework based on functional analysis. 
We finally emphasize that although the Hamiltonian with a Dirac's delta function potential, 
Eq. (pT|) , loses its meaning in the limit of N — > oo in two-dimensional billiards, Eq.(|42|) is 
still valid for any v, Vg if z < En with an arbitrarily fixed N. 

V. QUANTUM SPECTRUM OF BILLIARDS WITH A SMALL BUT 

FINITE-SIZE SCATTERER 

In this section, we discuss the quantum spectra of the billiards with a small but finite-size 
scatterer. We are after the wave-chaotic spectra in the systems with a realistic finite-size 
scatterer inside. For this purpose, we attempt to describe the finite-size scatterer in terms 
of a Dirac's delta function potential with a suitably truncated basis. This also enables us to 
recognize the physical meaning of the formal coupling constant vg. Relying on the arguments 
in Sees. Ill and IV, we clarify the condition for the appearance of wave-chaotic spectra for 
billiards with a small scatterer. 

Suppose that a finite-size scatterer of area Q is placed at x = xq inside an integrable 
billiard of area S. We describe the interaction in terms of a potential with a constant 
strength on the domain of the scatterer; 



case, since a oc p av (z) oc z 



1//2 , the RHS of Eq. (|37|) converges. In this case, one can relate v 



with Vg through 




(43) 
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V{x) = { (44) 
0, f G S - SI, 

where we denote the domains of the scatterer and the outer billiard by the same symbols as 
the areas. The matrix elements of V(x) are given by 

(V) mn = V/ ip m (x)tp n (x)dx, {m,n = 1,2,3, •• •). (45) 

./fi 

If the area of the scatterer is far smaller than that of the outer billiard, the scatterer is 
expected to behave as pointlike at low energy because waves with long-wavelength cannot see 
the shape of the scatterer. More precisely, the matrix element (V) mn with E m , i? n <(2Mr2) _1 
is suitably approximated by Vfl(p m (xo)(p n (xo) . Furthermore, the off-diagonal matrix ele- 
ments with E m <C (2Mf2) _1 <C E n are small, because of the mismatch in the wave numbers. 
These indicate that the energy spectrum of the billiards with a finite-size scatterer of area Q 
can be reproduced by the Hamiltonian Eq.([H]) in terms of a Dirac's delta function potential 
with strength 

v(Q) = VVL (46) 

along with a restricted basis up to 

E m) ~ (2MQ)- 1 . (47) 

It is crucial that the zero-range approximation is justified as long as z <C -EWn); even if |V| 
is sufficiently large. The mixing among unperturbed eigenstates is limited to a small number 
of states even if v (or V) is sufficiently large. Note also that the matrix V has a band-like 
structure, the width of which is independent of V, but is essentially determined by Q. 

We can now describe the low-energy spectra of the billiards with a finite-size scatterer 
inside within the framework in Sec. II. It is realized from Eq.(^) that, as long as z <C Enoi), 
the eigenvalues in case of finite-size scatterer can be calculated by Eq.fllSD with 



v e 1 ~ v(Sl) 1 + a ln(E N(Q) - z) 

~v{n)- l + a\nE m) . (48) 
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The meaning of the renormalized coupling constant vg becomes clear; A quantum billiard 
with a pointlike scatterer with strength vg can be obtained by taking the limit Q — > along 
with the limit 

" (Si) = B ,- + aL(2Affl) ~* "°- (49) 
or equivalently, with the limit 

1 

V = v^tt + attln(2Mtt) ~* ^ 



Any two sets of (V,fi), say (Vi,fli) and (V2,fi 2 ), which satisfy Eq. (|50|) with fixed v$, de- 
scribe the same low-energy dynamics as long as z <C £7/v(n)> where Q = max{f2 1 ,fi 2 }- The 
logarithmic term in the denominator of Eq. fl50|) expresses the reason why we cannot define a 
pointlike scatterer in terms of a Dirac's delta function in two dimension; The limit of Q — ► 
along with keeping VQ constant induces a potential which is too strong to define a quantum 
mechanical Hamiltonian for a pointlike scatterer. It is noteworthy that in the small-size 
limit, the strength of the scatterer is always negative in the sense that V — > — oo as Q — > 0. 
This is consistent with the fact that a single eigenstate with an eigenvalue smaller than E\ 
exists for any vg. We emphasize that this is not due to a specific formulation discussed here, 
but is required from the self-adjointness of the Hamiltonian. It should be noted, however, 
that as long as the obstacle has a substantial size, we can deduce all the low-energy physics 
from the present formulation by identifying v(Q) with vg through Eq. (f48D , even in the case 
of positive V. 

Inserting Eq.fl48|) to Eq. (|29|) , we obtain the condition for the appearance of the wave- 
chaotic energy spectrum; 



v(tt)~ l - aim 



En{q) — z 

Keeping Eq.(|23|) in mind, we can rewrite Eq.([H|) as 



<A~a. (51) 



'((vur 1 



In 



Pav J E N{n) - Z 



<1. (52) 
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Because the logarithmic function in Eq.(|51|) is negative for z < £7jv(n)/2, we obtain an 
interesting spectral property for lowest eigenstates; The wave chaos is most visible when 
the potential is attractive. As the energy increases up to z ~ E^iq\/2, the wave-chaotic 
region is expected to shift to — > according to Eq.([H]). However, the zero-range 
approximation tends to lose its validity at high energy, particularly in attractive cases as 
shown later in the numerical calculations. 

The previous numerical observations of Sinai's billiard give a supplementary evidence 
for the above finding. In Ref. [26|, Berry has exhibited the energy- level diagram of Sinai's 



billiard as a function of the radius of a circular obstacle in the billiard. It reveals that the 
integrability is restored in the low energy region as the radius decreases. This can be easily 
understood in the present scope. When the potential is repulsive, the wave chaos never 
appears even at low energy even if the repulsion is sufficiently strong. In spite of the infinite 
height of the potential in Sinai's billiard, the limit of the radius being zero can be taken 



without any difficulty because only the 'desymmetrized states' are considered in Ref. p6 |. 

In order to confirm the validity of Eq. Q5"T|) , we examine the quantum spectra of a rect- 
angular billiard with a small rectangular scatterer inside, which is analogous to the system 
discussed in Ref. || . The unperturbed eigenvalues and the corresponding normalized eigen- 
functions of the empty billiard are given by 

B ™ = 2M«-' + V >' (M) 

and 



/ 4 innx . niry ,_ AS 
<Pmn {x, y ) = J yy sin — — sin —— , (54) 

with m, n — 1, 2, 3, • • ■, respectively. The side lengths of the (outer) rectangle are denoted by 
l x and l y . Suppose that a rectangular scatterer is placed inside the outer rectangle such that 
the sides of the inner and outer rectangles are parallel to each other. We denote the side 
lengths and the position of the center of the scatterer by 5l x , 5l y and (xo,yo), respectively. 
The matrix elements of the potential of Eq.(|44]) are given by 



C^O m i ,rt i :rti'2 .n-i Vu mi rn2 {J-xi Slxy -£o)^ni,n2 (Jy > Sly j l/o) j (55) 

where 



2 r 1 (m— nWxn • (m— n)irSl 

-\ — — cos * H — 11 sm i — ^ — 

7r ^ m—n i 21 

1 (m+n)7ra;o • (m+n)7r<5£ 



m+n 



cos sin »}, m ^ n, (56) 



f - ^ cos sin 2^, m = n. 



In the actual numerical calculations, we set M = 2ir, l x = 7r/3 = 1.04719 and l y = 3/tt = 

0. 954929, which leads to p av = 1 and a ~ 1. We also assume Sl x = 3.53830 x 10~ 2 , 
Sly = 3.14023 x 10~ 2 and x = (x ,y ) = (0.622482,0.275835). In this choice, the area of 
the scatterer, Q, is 1/900 of the outer rectangle. 

Fig. 1(a) shows several lowest eigenvalues of this system for various values of the in- 
verse physical strength. The circles indicate the exact eigenvalues obtained by diag- 
onalizing the Hamiltonian matrix along with Eq. (|55"D . The strength of the scatterer, 
V, is determined by Eq. fl4"6|) . The value of v^ 1 ranges from —10 to 10 at intervals of 

1. The solid curves are the eigenvalues obtained by the approximation in terms of a 
Dirac's delta function potential. In the latter case, we truncate the unperturbed basis 
at E N( si x ,8i y ) = {(n/25l x ) 2 + (n/2Sl y ) 2 }/(2M) = 355.949, irrespective to the value of v. This 
corresponds to making a truncation at the momenta k x ~ tt/(2SI x ) and k y ~ tt/ (2Sl y ). ^From 
Fig. 1(a), we observe that the approximation of a zero-range potential with a suitably trun- 
cated basis is quite satisfactory, even if v is quite large. Since we have v ~ ((V) m ,n;m,n) / Pav 
with a ~ 1, = (y = 1000) indeed corresponds to a strong repulsion in the sense 
that the matrix elements of V are far larger than E^^i x ,si v )- Exceptions are the cases of 
a strong attractive potential. In this case, there is a considerable discrepancy between the 
exact and approximated eigenvalues. Fig. 1(b) shows the eigenvalues for large \v\. When the 
scatterer is attractive, there exist, in general, several eigenstates with negative eigenvalues; 
For example, three eigenstates with eigenvalues —1098, —261 and —151 for v^ 1 = —0.5, 
and four eigenstates with eigenvalues —2182, —1131, —944 and —12 for v^ 1 = —0.3. (These 
values, particularly for shallow eigenstates are not fully convergent despite a huge matrix 

19 



diagonalization with 12000 unperturbed states.) The eigenvalue of negative eigenstates is 
sensitive to the value of v. This causes a rapid change of the lowest positive eigenvalues as 
v varies. 

As shown in Fig. 1(a), the inflection points of the solid curves in the low energy region 
appear in a weakly attractive region. The prediction of Eq. ( pTj ) shown by a strip between two 
broken lines in Fig. 1(a) reproduces the position of the inflection points. In Fig. 2, we show 
the contour plot of the exact eigenfunction for several eigenstates. It can be seen that the 
mixture of the unperturbed eigenfunctions indeed occurs around the inflections points, while 
the wavefunction does not differ substantially from one of the unperturbed eigenfunctions in 
the region far from the inflection points. We stress that, even if \v\ is sufficiently large, the 
effects of high momentum components are visible mainly in the vicinity of the scatterer. The 
global behavior of eigenfunctions gradually changes along the solid curves as v _1 varies. It 
is expected that, as the scatterer shrinks (Ejsr(si Xi si y ) becomes large), the wave-chaotic region 
spreads over to higher energy. Moreover, the effective strength of the scatterer becomes 
weaker with smaller value for \v\ for a fixed energy in z <C En^ sl,), as indicated in Eq. (pl|) . 

Finally, we examine the validity of the zero-range approximation at higher energy. Figs. 3- 
5 show the eigenvalues in three cases; z ~ 0.10E N ^i x ,si y ) in Fig.3, z ~ 0.30Ejsr(si x ,5i v ) m 
Fig. 4 and z ~ 0.55^^(5;^,^^^) in Fig.5, respectively. It is observed that the accuracy of the 
zero-range approximation depends on the sign of the potential. For attractive cases, the 
zero-range approximation fails at somewhat low energy, as pointed out before. There is 
a considerable discrepancy even for v ~ — 1 at z ~ 0. 10 E^^i x ,si y )- On the other hand, 
for repulsive cases, the approximated eigenvalues are in fairly good agreement with the 
exact eigenvalues up to z ~ 0. 30 E N ^i x ,si y )- As the energy increases further, the zero-range 
approximation starts to fail except for in case of a very weak potential. It is easily recognized 
that the spectral properties at high energy is subject to the geometry of the scatterer, since 
a quantum particle moving in the billiard tends to see the shape of the scatterer as the 
energy increases. 
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VI. CONCLUSION 



We have examined the quantum spectra in two-dimensional billiards with a small-size 
scatterer inside. First, we have clarified the spectral properties of pseudointegrable billiards 
with a pointlike scatterer from a general point of view. The strength of a pointlike scatterer 
is specified by a renormalized coupling constant vg, formally defined within the formulation 
based on functional analysis. Although the coupling constant does not have a direct rela- 
tion to physical observables, it can be related to the physical coupling constant v defined 
as a strength of a Dirac's delta function potential together with a truncated basis. In two- 
dimensional billiard problems, the zero-range interaction cannot be rigorously described in 
terms of a Dirac's delta function. Nevertheless, v can be related to vg as long as the unper- 
turbed basis is truncated at a finite number. From a physical perspective, the truncation 
of the basis is a consequence of finiteness of the size of the scatterer, as the uncertainty 
principle implies. Applying the findings in pseudointegrable billiards to the cases of a small 
but finite-size scatterer, we have shown that the signature of wave chaos in lowest eigen- 
states are observed most prominently when the scatterer is weakly attractive. The numerical 
experiments of a rectangular billiard with a small rectangular scatterer inside corroborate 
these arguments. 
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FIGURES 

FIG. 1. Several lowest eigenvalues are shown for various values of the inverse physical strength 
of the scatterer, v~ l . The circles indicate the exact eigenvalues obtained by diagonalizing the 
Hamiltonian matrix along with Eq.(55). In this case, V = v/Q with Q = 1/900. While the value 
of v -1 ranges from —10 to 10 at intervals of 1 in (a), it ranges from —1 to —0.3 and from to 1 
at intervals of 0.1 in (b). The case of v~ l = is calculated by taking v = 1000. The solid curves 
are the approximated eigenvalues in terms of a Dirac's delta function potential with strength v, 
where the unperturbed eigenfunctions are truncated at E N ( &lxtSly -} = {(ir/25l x ) 2 + (ir/25l y ) 2 }/(2M), 
irrespective to the value of v. Unperturbed energies are indicated by vertical lines on the z axis. 
The strip between two broken lines in (a) is a prediction of Eq.(51) with a = 1, on which the 
signature of wave chaos appears. 

FIG. 2. The contour plot of the wavefunction for the eigenstates indicated by a-e in Fig.l is 
shown in (a)-(e), respectively. The value of for a-e is (5.0,4.54), (0.1,4.93), (-0.3,4.43), 

(—4.0, 5.63) and (—10.0, 6.13), respectively. The location of the scatterer is indicated by a rectangle 
in the billiard. 

FIG. 3. The eigenvalues around z ~ 0.10E N ( S i xSly y The indications are the same as in Fig.l. 
FIG. 4. The eigenvalues around z ~ 0.30E N ^i x ,si y )- The indications are the same as in Fig.l. 
FIG. 5. The eigenvalues around z ~ 0.55E N ^i Xt si y )- The indications are the same as in Fig.l. 
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